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Abstract. The recent COSY- 1 1 collaboration measurement of the two-proton correlation function 
in the pp — > ppr\ reaction, reported at this meeting [1], arouse some interest in a simple theoretical 
description of the correlation function. In these notes we present a pedagogical introduction to the 
practical methods that can be used for calculating the correlation function. 

Keywords: particle correlations, nucleon-nucleon interaction 
PACS: 25.75.Gz, 13.75.Cs 

We are going here to deal with low-energy phenomena and in order to avoid unneces- 
sary complications our approach will be non-relativistic. We wish to develop a practical 
scheme for calculating the two-particle correlation function C[k), with k denoting their 
relative momentum, choosing as the departure point the familiar formula 

C(k) = |D(r)|4>,(r)| 2 d 3 r, (1) 

(referred to in the literature as the Koonin-Pratt model [2]), expressing C(k) as an overlap 
of two distributions. The first distribution, D(r), is the effective source density function 
that is usually for convenience assumed to be a Gaussian 

D(r) = (4nd 2 )- 3/2 z~ {r/2d)2 (2) 

with a single parameter d reflecting the size of the source. The second factor entering the 
overlap integral (1), is a probability density involving the square of the wave function 
<&k( r ) describing two-particle system in the continuum. 

Actually, formula (1) is only a static approximation to the correlation function derived 
under the assumption that the final-state interaction between the two detected particles 
dominates, while all other interactions are negligible. Furthermore, it is assumed that the 
correlation functions are determined by the two-body densities of states (corrected for 
their mutual interactions) and that the single particle phase space distribution function of 
the emitted particle varies slowly with momentum. Admittedly, the question concerning 
the validity of these assumptions is far from settled but since they result in a manageable 
calculational scheme it is worthwhile to examine in some detail its consequences. 

Adopting hereafter units in which fi = c = 1, for non-interacting particles ^(r) takes 
the form of a plane wave 

* 4 (r) =e* r 

and in this case the correlation function is equal to unity, which means no correlation. 
However, when the particles are non-interacting identical bosons(fermions), the plane 



wave needs to be properly symmetrized(antisymmetrized), viz. 



*f (r) = 



q 



e ik r ±e ik r 
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and the interference term yields a non- vanishing contribution to the correlation function 
in excess of unity. For a Gaussian source, direct integration gives 

C ± (k) = \±e- 4d2k2 (3) 

and in the symmetric case C+(k) peaks at threshold with C+(0) = 2 going eventually 
to unity at large k. By contrast, in the antisymmetric case C (k) vanishes at threshold 
reaching unity from below for large k. The rate at which the correlation function ap- 
proaches unity is controlled by d, the only parameter in the game. Actually, C±(k) could 
be viewed as a function depending on two parameters k and d but for free propaga- 
tion, owing to dimensional scaling, the correlation function can only depend upon their 
product dk bearing a universal character. For interacting particles, as soon as additional 
parameters become available, further dimensionless quantities can be formed and dk is 
no longer the only possible combination. 

When the two-particle system has additional degrees of freedom the even and odd 
components of the wave function may both enter the correlation function. This happens 
in the two-proton case where the isospin part of the wave function is necessarily sym- 
metric and therefore Pauli principle admits in the spin- singlet states only even £ whereas 
in spin-triplet states, respectively, odd I. In the simplest case of two non-interacting pro- 
tons the resulting spin weights are purely statistical, and we get 



C pp (k) = \c + {k) + 3 -C-(k) = l- \^ ( X) 



with the intercept C pp (0) = j. 

Apart from correlations associated with permutation symmetry there would be also 
dynamical correlations induced by the two-particle interaction. To get some insight into 
the nature of dynamical correlations let us consider the case of two different particles 
whose propagation is not free but distorted by their mutual interaction. For simplicity 
we assume that the interaction may be represented by a spherically symmetric potential 
V(r). The wave function may be expanded in partial waves 



<D*(r) = £(2£+l)iVM/(^(M (5) 

where Pg(z) denotes Legendre polynomial and the function ug(k,r) is the solution of 
the appropriate radial Schrodinger equation and for r — > °° satisfies the outgoing wave 
boundary condition 



ug(k, r) ~ sin(£r — + kfg(k) exp [i(kr — j£n)] 



(6) 



where fe(k) is the partial wave scattering amplitude. Inserting (5) in (1), after trivial 
angular integration, the correlation function takes the form of a series 



Ait 00 r 00 
C W = T2E( 2 ^+ 1 )/ D(r)\u e (k,r)\ 2 dr (7) 

whose convergence rate depends upon the case considered. For a short ranged potential 
the functions \ue(k, r)\ 2 are pushed outward by the centrifugal barrier and their overlap 
with the density rapidly decreases with increasing £ and in consequence the series (7) 
also rapidly converges. The same reasoning repeated for the symmetric wave functions, 
yields 

C +( k ) = TI E ( 2£+1 ) / D(r) \u e (k,r)\ 2 dr (8) 
K V=even Jo 



and for C (k) analogous formula holds but the summation would be over odd £. 

As an important example we consider the Coulomb interaction V c (r) = a/r where a 
is the fine structure constant times a product of the charge numbers. The radial functions 
are well known 

u e (k,r)=F e (ri,p) (9) 

with p = kr, r\ = (Xfl/k where /1 is the reduced mass and the regular Coulomb 
functions Fi(r\,p) are defined in Abramowitz and Stegun [3] (for numerical methods 
of computing them cf. [4] and references therein). If the two protons experienced only 
Coulomb interaction, the pp correlation function respecting Pauli principle would be 



7.TZ \ f°" 

C Pp( k ) = Tl{ E (2^+1)/ D(r)F e (l 1 ,kr) 2 dr + 

K [e=even J0 

f'OO 

+ 3 £ (2£+l) / D(r)F e (riM) 2 dr 



(10) 



where the integrals require numerical treatment. 

The last step is the inclusion of strong interaction which, in general, leads to compli- 
cations stemming from the fact that the nuclear forces have non-central components like 
the tensor interaction or spin-orbit force, and the orbital momentum is not a good quan- 
tum number. However, for low energies (k < 100 MeV /c, say) the dominant contribution 
comes from the s-wave and the interaction in higher partial waves can be neglected. In 
this case £ is still a good quantum number and including both, the Coulomb and the 
nuclear s-wave interaction, the pp correlation function respecting Pauli principle, reads 

C Pp( k ) = 2 £;\ rD(r)\u (k,r)\ 2 dr+ £ (2£+l) [°° D(r)F e (r,,kr) 2 dr \ + 

K { J0 £=even>0 J ° J 
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-j £ (2£+l) f °D(r)F t (ri,kr) 2 dr 

K i) 1,, JO 



where uo(k, r) is the solution of the wave equation for i = involving both, the Coulomb 
and the strong interaction potential. The above formula constitutes the basis for all our 
computations. 

To get a feeling how (11) works in practice it is instructive to provide some example. 
Perhaps the simplest simulation of the nuclear force delivers the delta-shell potential 
model when V(r) is taken in the form 

2flV(r) = -(s/R)8(r-R) (12) 

where s is a dimensionless strength parameter and R is the range parameter, assuming 
that the potential (12) is operative in s-wave only. The wave function uo(k,r) may be 
obtained by solving the appropriate Schrodinger equation but here we shall use for that 
purpose the equivalent Lippmann-Schwinger equation which incorporates the correct 
asymptotic boundary condition. More explicitly, we are presented with the integral 
equation 

f'OO 

u (k,r)=F (ri,kr)+2ll J g+(r,r')V(r')u (k,r')dr f (13) 
where g£ (r,r') is the outgoing wave Coulomb Green's function 

g+(r,r') = -(l/k)F (r],kr<) [Go(TUr>)+iF (77>>)] (14) 

with G((r\,kr) denoting the irregular Coulomb wave function (cf. [3], [4]) and the 
other symbols are: r< = min(r,r / ); r> = max(r,r / ). Note, that acting with the operator 
d 2 /d^ + k 2 — Ifia/r on both sides of (13), one immediately recovers the Schrodinger 
equation. For the potential (12) it is a trivial matter to obtain the solution of (13) and the 
latter is given in an analytic form 

(kr)= U(k)F Q (r],kr) for r<R 

U0[ ,rJ \F Q {7]M)+h{k)k[G Q {7]M)+m{7]M)] for r > R 

where the amplitude A (k) is 

A ^ = 1 - (s/kR)F (ri,kR) [G (ri,kR) +iF (ri,kR)} (16) 

and the second constant fo(k) which is the appropriate pp scattering amplitude, is given 
as 

fo(k) = (s/k 2 R)F (ri,kR) 2 A(k). (17) 

Feeding (11) with uo(k,r) given in (15), and performing (numerically) the integrations 
we obtain the correlation function. For "realistic" nuclear potentials the procedure would 
be similar, the only difference being that uo(k : r) must be generated numerically from 
(13). 

In our computations we have tried a variety of NN potentials (Reid68, Reid93, Ni- 
jmegen93 and Argonne AVI 8) but the corresponding correlation functions were in all 
cases nearly indistinguishable. On the other hand, the delta-shell potential resulted in 
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FIGURE 1. (left) Two-proton correlation function versus k for different source radii calculated from 
the Reid soft core potential, (right) Partial wave (£ — 0) pp cross-section calculated from (19). The stars 
show the cross-section computed from the toy model formula (17) for s = 0.906 and R = 1 .84 fm. 



markedly different shape of the correlation function. The plots of the correlation func- 
tion calculated for Reid soft core potential are presented in Fig. 1 (left) for different radii. 
In qualitative terms the curves are similar: all show a depression close to threshold and 
a prominent peak at about 20 MeV/c. The same behavior has been observed experimen- 
tally [1]. In our model the position of the peak does not seem to be dependent upon d 
whilst its height is very sensitive to the source radius. The depression close to threshold 
results from a combined effect of Fermi statistics and Coulomb repulsion which both try 
to keep the two protons apart causing the correlation to be negative at small k. The peak 
at 20 MeV/c has a dynamic origin, manifesting strong nuclear attraction experienced 
by the protons. The height of the peak decreases rapidly for increasing d assuming the 
largest possible value for d = 0. In the latter case the source shrinks to a point and the 
density (2) can be replaced by a delta function D(r) = 8(r). The correlation function for 
a point-like source can be immediately obtained, and reads 

taC^t) = \ ^ff- = ^ «<*,0)I J , (18) 

where prime denotes derivative with respect to r. The expression on the right hand side of 
(18) is nothing else but the enhancement factor associated with pp final-state interaction, 
used e.g. to approximate the pp — > ppf] cross section, and which is also known to have 
a prominent peak at 20 MeV/c. 

To understand better the origin of the peak at 20 MeV/c let us consider first a 
simplified situation with Coulomb interaction switched off. Close to threshold the NN 
t = scattering amplitude may be expressed in terms of the scattering length a as 
fo(k) = (— 1 /a — ifc) ~ 1 and in the complex k-plane thi s amplitude has a pole at £: p = i/a . 
For a > the pole is on the physical sheet (Zk p > 0) and we have a bound state. This is 
the deuteron case. However, when a < 0, the pole is on the non-physical sheet (3k p < 0) 
and we are left with a virtual state. For the pp pair this would be the particle unstable 



2 He state. In either case the NN cross section shows a pronounced peak at threshold as 
seen in Fig. 1 (right) (dashed curve). Putting back the Coulomb interaction and including 
the effective range (ro) term, a model independent expression for the low-energy pp 
scattering amplitude, takes the form 

f (]A ~ C o(^) ng\ 

-\/a + \rok 2 -2kT)h(T))-ikCl{T)) 

where Cq ( r\ ) = 2nr\ / [exp (2nr\ ) - 1 ] ) is the Gamow factor, h ( r\ ) = 9t y( 1 + ir\ ) - log | r\ \ 
and \y denotes the digamma function [3]. The corresponding partial wave cross-section, 
i.e. \ fo(k) | 2 for a=-7.78 fm and ro=2.72 fm is presented in Fig.l(right) by the continuous 
curve and a similar result would be obtained from the toy model formula (17) (depicted 
by stars in Fig. 1 (right)). Owing to the Coulomb corrections, the peak in the cross section 
is shifted off-threshold by about 20 Mev/c and its height is depressed to one third 
of its original value. This is now the very same peak that appears in the correlation 
function and can be regarded as an artifact of the 2 He state. It is not a resonance though 
because the real part of the denominator in (19) equal — 1 /a + \rok 2 — 2kr\h{r\) does not 
vanish, showing instead only a minimum whose position is identical with the position 
of the peak in Fig. 1 (right). Actually, given the low-energy pp scattering parameters, the 
position of the peak could have been predicted by equating to zero the derivative of the 
real part of the denominator in (19) and by solving the resulting equation for k. 
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